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ANISOTROPIES IN ULTRAHIGH ENERGY COSMIC RAYS* 

JOHN SWAIN 

Department of Physics, Northeastern University, Boston, MA 02115 
e-mail: john.swain@cem.ch 

The present status of anisotropy studies for the highest energy cosmic rays is presented 
including the first full sky survey. Directions and prospects for the future are also 
discussed in light of new statistical methods and the last quantities of data expected in 
the near future from the Pierre Auger Observatory. 

1. Introduction 

Ultrahigh ener^ cosmic rays (UHECR's) are among the most enigmatic phenomena 
in the universe.*^ In the mid-60's Greisen, Zatsepin, and Kuzmin (GZK) pointed out 
that ultra high energy protons interact with the all-pervading cosmic microwave 
background via photopion production.El Strictly speaking, protons with energies 
> 10^^'^ eV have a mean interaction length < 6 Mpc and an inelasticity of about 20% 
per interaction. Consequently, the popular astronomical picture, namely proton 
"bottom-up" acceleration in extragalactic objects, predicts a sharp suppression of 
the cosmic ray intensity somewhat beyond 10^^'^ eV. For heavy nuclei, the giant 
dipole resonance can be excited at similar total energies and hence iron nuclei do 
not survive fragmentation over comparable path lengths.l^ 

The existence of cosmic rays with energies exceeding 10^° eV has been observed 
by the Volcano Ranch, ^ the Haverah Park,^ the Sydney University Giant Airshower 
Recorder (SUGAR), ^ the Akeno Giant Air Shower Array (AGASA),'^ and the Fly's 
EyJ^experiments. Because of the GZK cutoff, these cosmic rays should be produced 
in nearby active astronomical objects. Such high energy "stars" have been searched 
for in the arrival direction of these events but no clear candidates were found. 

Of course, there are ways to avoid the distance restriction imposed by the GZK 
effect. For instance, there could be a "top-down" mechanism where (charged and/or 
neutral) supermassive A-particles are produced at extreme energies. Sources of 
these exotic particles could be topological defects left over from early universe phase 
transitions associated with the spontaneous symmetry breaking that underlies uni- 
fied models of high energy interactions,!^! or else some long-lived metastable super- 
heavy relic particles produced through vacuum fluctuations during the inflationary 
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stage of the universe.^l From time to time, the energy stored in a single X can 
be released in the form of massive quanta that typically produce jets of hadrons 
well above the highest observed energies. However, it is noteworthy that there can 
be a problem with some top-down interpretations. Specifically, the X-particle cas- 
cades may produce a rather large flux of energeti c ph otons and neutrinos, possibly 
in excess of the upper limits already established.l^ The lack of plausible nearby 
astrophysical sources has also encouraged the idea of positing undiscovered neutral 
hadrons,^! as well as mechanisms which are able to break the GZK barrier Even 
though sufficiently heavy particles would avoid photopion production (the threshold 
energy varies as the square of the mass of the first resonant state), the existence of 
these particles now appears to be excluded by laboratory experiments."'^'^ The only 
standard model particle that can reach our galaxy from high rcdshift sources with- 
out significance loss of energy is the neutrino. However, the expected event rate for 
early development of a neutrino shower is less than that of an electromagnetic or 
hadronic interaction by several orders of mag nitu de, even if ones takes into account 
non-standard graviton mediated interactions.^^ 

The distribution of arrival directions is perhaps the most helpful observable in 
yielding clues about cosmic ray origin. This may come either from clustering on 
a small angular scale that identifies discrete sources,^ or else as a large-scale ce- 
lestial pattern that characterizes a particular class of potential sources.^l Along 
these lines, data observed by the AGASA,CHlthe SUGAR,IEland the Fly's EyeEOl 
experiments show an excess flux of cosmic rays from a direction near the Galactic 
center up to about 10^^'^ eV, but there is evidence for Galactic plane avoidance 
above that energy. Such an effect can be easily explained if cosmic rays are mostly 
protons and nuclei, because their magnetic rigidity increases with energy and so 
one expects the angular width of the Galactic plane as seen in protons and nuclei 
would shrink slowly with rising energy. Moreover, the events yielding the observed 
anisotropy are concentrated in a limited energy range. This is very suggestive of 
neutrons as candidate primaries, because the directional signal requires relatively- 
stable neutral primaries, and time-dilated neutrons can reach Earth from typical 
Galactic distances when the neutron energy exceeds 10^^ eV. Arguably, if the Galac- 
tic messengers are neutrons, then those with energies below 10^^ eV will decay in 
flight, providing a flux of cosmic antineutrinos above a TeV which is observable at 
a kilometer-scale neutrino observatory!^ A measurement of this flux can serve to 
identify the first extraterrestrial point source of TeV antineutrinos. 

All in all, the data around 10^^'^ eV suggests that a new population of cosmic 
rays with extragalactic origin begins to dominate the more steeply falling Galactic 
population. In the extragalactic sway, the evidence for anisotropy patterns is sug- 
gestive but statistically very weak. On the one hand, correlations with the local 
structure of galaxies have been reported,l23 suggesting that all cosmic rays with en- 
ergies > 10^*'^ eV are emitted by nearby astrophysical sources. On the other hand, 
possible correlations with high redshift astrophysical objects, which most likely indi- 
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cate new physics, are also under debate!^ Clearly, a positive identification of ultra 
high energy cosmic ray sources from the distribution of arrival directions requires 
a careful study of such distribution over the full celestial sphere. The latter is the 
main inspiration for this talk, which is organized as follows. In order to set the 
stage for the full-sky coverage discussion, in Sec. 2 I will present an overview of 
available statistics and existing anisotropy studies. In particular, I will concentrate 
on the SUGAR and AGASA experiments and discuss in some detail the exposure 
of these ground arrays. Ne xt, in Sec. 3, I will review the main properties of the 
angular power spectrurJ^ and study cosmic ra y an isotropics from an expansion 
on spherical harmonics for modes out to ^ = 5.^' In Sec. 4 a numerical likeli- 
hood approach to the determination of cosmic ray anisotropies is presented. This 
method offers many advantages over other approaches: It allows a wide range of 
statistically meaningful hypotheses to be compared even when full sky coverage is 
unavailable, can be readily extended in order to include measurement errors, and 
makes maximum unbiased use of all available information. Finally, I will summarize 
and present the conclusions in Section 6. 

2. Experimental data sets 

The SUGAR array was operated from January 1968 to February 1979 in New South 
Wales (Austraha) at a latitude of 30.5° South and longitude 149°38' East.^^ The 
array consisted of 47 independent stations on a rectangular grid covering an area 
S" « 70 km^. The primary energy was determined from the total number of muons, 
Nfj^, traversing the detector at the measured zenith angle 9. The total aperture for 
incident zenith angles between 9i and 02 was found to be 



Here, p(iV^, 0) is the probability that a shower falling within the physical area was 
detected, S'cos6' is the projected surface of the array in the shower plane, and d^l 
is the acceptance solid angle. The SUGAR Collaboration reports^^ a reasonable 
accuracy in assessing the shower parameters up to 6* = 73°. The estimated angular 
uncertainty for showers that triggered 5 or more stations is reported as 3° sec 
However, the majority of events were only viewed by 3 or 4 stations, and for these the 
resolution appears to be as poor as 10°. -2^ Of particular interest for this analysis, 
p{Nf, > 10^6I < 55°) w 0.85,Eni yielding a total aperture A w 125 km^ sr. This 
provides an exposure reasonably matched to that of AGASA, which is described 
next. 

The AGASA experiment occupies farm land near the village of Akeno (Japan) at 
a longitude of 138°30' East and latitude 35°30' North.lSIl The array, which consists 
of 111 surface detectors deployed over an area of about 100 km^, has been running 
since 1990. About 95% of the surface detectors were operational from March to 
December 1991, and the array has been fully operational since then. A prototype 
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detector operated from 1984 to 1990 and has been part of AGASA since 1990.1211 
The aperture for events with primary zen ith angle 0° < 6 < 45° and energies beyond 
-|^gi9.25 gY jg found to be A « 125 km^ sr."^^' The angular resolution for these events 
is 1.6°. El 




A detector at latitude oq that has continuous operation with constant exposure 
in right ascension and is fully effici ent for 9 < ^max has relative exposure with the 
following dependence on declinatiorl^ 

Ci;((5) cx (cosao cos 6 sin amax + Q^max sinao sin (5) , (2) 

where amax, the local hour angle at which the zenith angle becomes equal to ^max, 
is given by 

if ^ > 1 

TT if C < -1 (3) 

cos~^ ^ otherwise 

with 

_ cos ^niax " sin oq sin S 



cos ao cos S 



(4) 



The resulting declination dependence for SUGAR and AGASA together with the 
combined aperture is given in Fig. ^ 

As one can readily see in Fig. ^ the combined aperture of SUGAR and AGASA 
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arrays is nearly uniform over the entire sky. The expected event rate is found to be 

1 1 



El El 



(5) 



where {E^J{E)) w lO^'*-^ eV^ m ^ s ^ sr ^ stands for the observed ultra high 
energy cosmic ray flux, which has a cutoff at E2 = 10 Thus, in approxi- 

mately 10 yr of running each of these experiments should collect ~ 50 events above 
El = 10^^-^ eV, arriving with a zenith angle < 6'max- Here, 6'niax = 45° for AGASA 
and ^niax — 55° for SUGAR. Our sub-sample for the full-sky anisotropy search 
consists of the 50 events detected by AGASA from May 1990 to May 2000,'^ and 
the 49 events detected by SUGAR with 6 < 55°. Note that we consider the fuU 
data sample for the 11 yr lifetime of SUGAR (in contrast to the 10 yr data sample 
from AGASA). This roughly compensates for the time variation of the sensitive 
area of the experiment as detectors were deployed or inactivated for maintenance. 
The arrival directions of the 99 events are plotted in Fig. |21 (equatorial coordinates 
B.1950). 




Figure 2. Arrival direction of the 99 events observed above 10^^'^ eV by the SUGAR (6* < 55°) 
and the AGASA {8 < 45°) experiments (equatorial coordinates B.1950). 



3. Correlations and Power Spectrum 

We begin this section with a general introduction to the calculation of the angular 
power spectrum and the determination of the expected size of intensity fluctuations. 
The technique is then applied to the AGASA and SUGAR data in order to check 
for fluctuations beyond those expected from an isotropic distribution. 
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Let us start by defining the directional phase space of the angular distribution 
of cosmic ray events in equatorial coordinates, {a, S). (i) The direction of the event 
is described by a unit vector 

n = sin (5 (i cos a + j sin a) + k cos S ; (6) 

(ii) The solid angle is given by 

(fn = sin 5 dSda ; (7) 

(iii) The delta function for the solid angle is defined as 

(S(n, n') = (5(cos(5-cos(5') ^ S{a ~ a' + 2Trm) , (8) 

m— — oo 

so that, as usual, 

/(n)= J S(n,n') f{n') d'n' ; (9) 

(iv) The probability distribution P(n)(i^n of events can be employed for the purpose 
of computing the averages 

7= I f{n) P{n) d'n; (10) 

Finally, (v) for a sequence of N different cosmic ray events (rii, . . . ,njv) one may 
assume an independent distributions for each event, i.e. 

N N 

P^(ni,...,njv)n^'n, = J]{P(n,)d'n,} . (11) 

i i 

For a sequence of events (ni, . . . , un) let us describe the angular intensity as 
the random variable 

1 ^ 

From Eqs. (^TJ and it follows that 

/R = J ... J Iin)PNini,...,nN)Y[d^n, 

= P{n)- (13) 
The two point correlation function G{n, n') = /(n)/(n') is defined via 

G(n, n') = / • ■ • / ^(n) Hn') Pjv(ni, . . . , njv) J] '^""^ 



= ls{n,n')P{n)+fl-^] P(n) P(n') 



(14) 
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The "power spectrum" of the correlation function is determined by the eigen- 
value equation 

J G(n, n') VaK) d^n' = A VA(n). (15) 
In this regard it is useful to introduce Dirac notation to indicate the inner product 

(^\^) ^ J ^l;*{n) V'(n) d^n . (16) 

With this in mind, Eq. H15|l reads 

G IVa) = A \,Px} . (17) 
In the limit of a large number of events iV — )■ oo, 

lim G(n,n') = Goo(n,n')-P(n)P(n'), (18) 

N—*oo 

or equivalently, 

Goo = IP) {P\ ■ (19) 

In such a limit, fluctuations can be neglected and we find only two possible values 
in the spectrum: (i) There is a non-degenerate non-zero eigenvalue 

Goo \P) = Aoo \P) , (20) 

with 

Aoo = {P\P) = J P^in)d^n. (21) 

(ii) For every state |/) orthogonal to \P) with mean value / = {P\f) — 0, there 
exists a zero eigenvalue in the power spectrum 

Goo I/) = \P) (P\f) - 11 P/d'n| \P) = / |P) = . (22) 

Let us now turn to consider the effects of finite N. Defining the fluctuations in the 
intensity by 



A/(n) = /(n) - /(n) = /(n) - P(n), (23) 
the two point correlation function can be re-written as 



G(n,n') = /(n)/(n') = /(n) I{n') + A/(n)A/(n') 

= Goo(n,n') + A/(n)A/(n'), (24) 

with 



A/(n) A/(n') - ^ [ '5(n, n')P(n) - P(n)P(n') ] , (25) 

where Eq. (|14|l has been invoked. Putting all this together, some general results 
follow: (i) For the N ^ oo case, there is only one state with a finite eigenvalue 
Aoo, while the rest of the power spectrum corresponds to A = 0. (ii) For finite N, 
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Eq. H25|l implies that the fluctuations are of order N~^^. The power spectrum for 
large N then has one eigenvalue of order unity and the rest of the eigenvalues are 
of order N^^. 

Now, for an isotropic distribution of n, 

Pin) = ^ , (26) 
and the two point correlation function reads, 

S(n, n') + - — -T7 1 

4ttN ^ ' ' (47r)2 V N 



GM = TZTr + 773^2 ( 1 - - ) ■ (27) 



The eigenvalue problem is solved by employing spherical harmonics 1^ 

J G(n, n') Y,min') d^n' = A,™ y,,„(n) , (28) 

where 

X _r(4^)-' if (£,m) = (0,0) 

^'"^ ~ \ (4^iV)-i if {£, m) ^ (0, 0) ■ ^ ' 

The eigenfunctions form a useful set for expansions of the intensity over the celestial 
sphere 

oo i 

/(n) = ^ ^ ai„iY£„i{n) . (30) 
e=o m=-e 

To incorporate the dependence on declination given in Eq. let us re-define the 
angular intensity 

1 ^ 1 

where ujj is the relative exposure at arrival direction nj and M is the sum of the 
weights ujj^. Since the eigenvalues of the Y^m expansion are uniquely defined 



J /(n) r,™(n) d^n , (32) 



the replacement of Eq. H31|l into Eq. H32(l leads to the explicit form of the coefficients 
for our set of arrival directions 

1 ^ 1 

^^"^ = Jf^ll ~ ^i"i{nj) . (33) 

j=l 3 

With the coeffients given in this way, one can plot the intensity of the cosmic ray 
sky using Eq. 1301 as seen in Fig. |2| The mean square fluctuations of the coefficients 
are determined by the power spectrum eigenvalues according to 



-,2 _ 



(34) 
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Right Ascension (deg.) 



Figure 3. Intensity of the cosmic ray sky in equatorial coordinates as seen by the SUGAR and 
AGASA arrays. 



Although fuU anisotropy information is encoded into the coefficients aim (tied to 
some specified coordinate system), the (coordinate independent) total power spec- 
trum of fluctuations 



Ci£) 



{2£+l) 



(35) 



provides a gross summary of the features present in the celestial distribution to- 
gether with the characteristic angular scale(s). Note that Eqs. H29|l and (|34|l imply 



1 



(2^+1) 



E 



(47r)"i if£ = 
(47riV)-i if £ ^ 



(36) 



The power in mode £ is sensitive to variation over angular scales of £~^ radians!^ 
Recalling that the estimated angular unc erta inty for some of the events in the 
SUGAR sample is possibly as poor as 10° 1^ we only look in this study for large 
scale patterns, going into the multipole expansion out to £ — 5. 

Our results at this juncture are summarized in Fig.0] The angular power spec- 
trum is consistent with that expected from a random distribution for all (analyzed) 
multipoles, though there is a small (2i7 ) ex cess in the data for £ — 3. The majority 
of this excess comes from SUGAR data.'i"' The decrease in error as £ increases may 
be understood as a consequence of the fact that contributions to mode £ arise from 
variations over an angular scale £~^. If one compares to the expectation for isotropy, 
structures characterized by a smaller angular scale, and hence larger £, can be ruled 
out with more significance than larger structures. 

To quantify the error, we study the fiuctuations in C{£) for £ > I. For simplicity, 
let us neglect the small effects of declination (viz., oJi = 1 Vi), and consider the 
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Figure 4. The angular power spectrum is indicated by the squares. The horizontal lines indicate 
the mean value, C{£) = (4:1tN)~^, expected for an isotropic distribution. The upper shaded band 
shows the la fluctuation around the mean value for N = 99. The band was obtained from 1000 
sets of Monte Carlo simulations of 99 events each, including small corrections for Wi. For £ = 3, 
where there is a small excess compared to the expectation for isotropy, C(3) = 2.16 X 10"'^ while 
the expectation from a random distribution is Cmc(3) = 9.5 X 10"*, with a variance of 5.0 X 10~^. 
The projected sensitivity for the Pierre Auger Observatory is also indicated on the plot by the 
lower shaded band. 



random variable 



C{£) ( Attn . ^ , 

Denoting by Pe(cos6) the Legendre polynomial of order £ and employing the addi- 
tion theorem for spherical harmonics, 



47r 



2i- 



- J2 YeMYe„,{n') - P,(n • n') , (38) 



Eqs. (EZIl, and (I2HIl imply that 



^^ = 1 + ^ E ^^n, -n,). (39) 



N 

l<i<j<N 
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Evidently, Xi — 1. Besides, 

4 

7V2 



^l^^ + 4^ E E Pdn,-n,)Piin,-n,). (40) 

l<i<j<N l<k<q<N 



Since different pairs in the sum on the right hand side of Eq. H40f) arc uncorrclated, 
it follows that 

^=1 + ]^ E ^Hn.-n,)2. (41) 

l<i<j<Ar 

There are {N{N — l)/2} equivalent pairs in Eq. 141(1 which implies 



X^^Xi +2(1--) Piim ■ n2)2 . (42) 



From Eq. H38|l we obtain 



P,(ni . n,r = (^^^ j E 
1 



m— — l 



2e + i' 

Plugging Eq. (gSl into Eq. 1321) leads to 



1 



or equivalently. 



yielding (for large N) 



which is the variance on Xf 



(43) 



hm ^ = for £>1, (46) 



4. Numerical Likelihood Analysis of Cosmic Ray Anisotropies 

The approach described above requires what might be described as "reasonably 
good statistics". That is to say, one needs to have essentially some reasonable 
acceptance for every element of solid angle in the sky. Problems of interpretation 
can easily arise if there are regions which are unobservable - blind spots - and in 
this section I briefly describe h ow t hese can be dealt with numerically. The work 
follows closely that presented inEHI 
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4.1. Spherical Harmonics 

It is convenient to define a set of real spherical harmonics {V'^,™,} for £ = 0,1,2, .. . 
and m = -£,-{£-!),. ..,0,. ..,{£- l),£hj 

^|m|pjm| ^^^g cos(m0) for m = —£, . . . , — 1, 

^l,m{0,4') = { fcOpO(cos6l) for m = 0, (47) 

kf^Pp{cos 6) sin(m(/)) for rr?, = 1, .... ^ 

which are orthonormal with respect to the usual measure sui{6)d9d(j) and integration 
over 6 from to tt and cp from to 2n, the are associated Legendre polynomials 
of the first kind, and the normalization constants kL, are: 



The natural measure of anisotropy for a spherical distribution is in terms of these 
spherical harmonics, as each £ labels, in a coordinate-independent fashion, just how 
much of each irreducible S0(3) representation is present. 

In a perfect world with infinite statistics and complete sky coverage there are 
now many possible approaches to estimating how much of each of these compo- 
nents is present in a distribution, or, better, what is the likelihood that a given 
function f{9, 0) with Fourier-Legendre expansion f{0, </>) = YleLo 'Y^ln--t '^i,m,''Pe,m 
representing the probability density of sources gives rise to the observed distribution 
g{9, (j)). The coefficients can be extracted from the usual integral 



ai,m= / / f{9,,p)^Pi,m{0,<P) sin 9d9d,p (49) 
Jo Jo 

Of course this gives no measure of what sort of error should be associated with 
the determined values of each coefficient. Alternative approaches are to fit for the 
coefficients by minimizing some x^-like quantity like 

/ ^ sm9d9d(j) (50) 

Jo o- 

with a a suitable measure of error, or to compute and maximize a corresponding 
likelihood that the hypothesized distribution parametrized by the ae^m gives rise 
to the the observed distribution g. Of course this is statistically unreasonable for 
a finite number of sources {i.e. to provide an infinite number of coefficients!). In 
general a decision must be made to truncate the expansion at some value of £ to 
make the sum finite, but now the general phenomenon of aliasing risks that, for 
example, a fit allowing for £ — 0^1 might give misleading results for observed data 
which is drawn from a purely £ = 2 distribution, say. 

A more serious problem is that should part of the sky be unobserved, there 
is now no way to calculate anything! This is not a trivial point. An attempt to 
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find anisotropies based only on observations in the Northern hemisphere with ze- 
roes inserted for the whole Southern hemisphere would be wildly in error if some 
simple extrapolation were made to the unobserved region of the sky - especially if 
there were something bright and as-yet undetected in the South! The real challenge 
is to say something statistically meaningful with the data that is actually avail- 
able. Clearly, observing only the Northern hemisphere and seeing a good degree 
of isotropy should increase one's net belief in overall isotropy of the full sky, while 
leaving open the possibility of a staggeringly bright or empty sky in the South. 
One approach sometimes advocated is to try to make some sort of new "special 
functions" which would be orthonormal over the observed region of the sky, but 
it's not clear that this has much physical meaning as it elevates a defect like "lack 
of acceptance" to a status comparable to "the S0{3) invariance of space embodied 
in the V'f ,m" • The following is a proposal for what seems to make good statistical 
sense, and is physically unbiased. 

4.2. A Likelihood Proposal to Handle Limited Acceptance 

Based on the above observations, the following proposal seems reasonable: keep the 
spherical harmonics as always with the (necess arily ) truncated Fourier-Legendre 
expansion and construct an unbinned likelihood^ function in which one clearly 
specifies which values oi£ are included in the sum. An unbinned likelihood function, 
as described below, automatically makes maximum use of all detected information, 
allows for measurement errors to be included easily, and is easy to implement nu- 
merically. Most importantly, however, the likelihood is not to be normalized as it 
stands. Rather, one should take the Bayesian approach to likelihood which says 
that likelihoods give us ways to update our prior experimental data or guesses in 
light of new information. This will mean that one can present results on various 
hypotheses about data without bias, and with the easy inclusion of other data from 
the same, or other experiments. 

To be concrete, we specialize here to the case where g is a sum of delta func- 
tions representing sources i of unit intensity at 9i , 4>i and later discuss how one can 
treat the case of these being at uncertain locations, or taking into account other 
properties such as intensity, energy, composition, etc.. These will appear as natural 
generalizations to the approach described. 

The (unnormalized!) likelihood L{{£}\9i, (f)i) that the measured 9i, (pi arise from 
f{0,(p) = ^m=-e'^i,"^'^^,"^ where the sum over £ is specified according to 
whatever hypothesis is being tested {i.e. just taking {£} = {0, 1} allows for uniform 
and dipole contributions and no others, while {£} — {2} would be pure quadrupole) 
is 

L{{£}\9,, <t>i) = (51) 

Jill nm=-£ dat,mf{9,(t)) 

where the integral in the denominator is over all the parameters that can vary and 
over the range in which the parameters are allowed to vary. An important caveat is 
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that for / to represent a sensible probability distribution it should never be negative, 
and this must be checked. Two approaches are possible: one is to restrict the domain 
over which coefficients range so that the function is strictly positive (this is not 
actually very difficult in practice since distributions are often nearly uniform with 
small fluctuations superimposed) or to take as a probability distribution function 
some positive function of / in place of / above. In astronomy [4] it is not uncommon 
to see exp(/). The choice ultimately represents the unavoidable presence of some 
(often hidden) assumption about what a sensible prior {i.e. in the absence of data) 
is for the likelihood - a point to which we return later. 

By construction then L{{£}\9i, (pi) is normalized so that its integral over all the 
parameters that can vary (the coeflacients ae^m included in the truncated Fourier- 
Legende expansion), and the total likelihood Ltot, which is a function of those 
ae,m is 

LTOTm) = ll^m\0u<l>i) (52) 

i 

This is a purely relative likelihood, and while no absolute normalization is possible 
(nor should it be!) if part of the sky is unobserved, it is now very useful for two types 
of calculation. If one has a prior expectation for the distribution of the (which 
might be that they are all a priori equally likely) then Ltot can be multiplied by 
this and the product treated as a normalized likelihood distribution for the o^,m 
themselves. This is, of course, potentially dangerous, but does allow one to see how 
the new data (the measured {di, (f>i)) should cause one to revise earlier beliefs. Such 
a likelihood can be maximized with respect to the a^,m aud if not enough data is 
present to test the hypothesis (for example, more coefficients to fit for than data 
points) the fit will respond by just not converging (that is to say, there will be flat 
directions in the likelihood as a function of the parameters meaning that one can't 
decide) - the beauty of this approach is that it is, by construction, correct. When 
results arc obtained, one automatically gets the values of the parameters, and their 
entire likelihood distributions from which errors (which need not even be Gaussian) 
can be extracted. One can even do things like fix some parameters to those given 
by a favoured theory and then repeat the fitting and then obtain likelihoods for the 
correctness of that theory. 

More objectively, one can compute relative likelihoods in which the prior drops 
out, so that it is reasonable to ask (even in the absence of full sky coverage!) 
what the relative likelihood Lrel is of pure dipole distribution to one admitting 
uniform, dipole and quadrupole components: Lrel = LTOT({o^ii]) '-"^'^ wants 
statements made to be about all energies over x, then one just uses the data points 
with energies over x. Similarly, data can be selected by composition and (even 
relative) anisotropies be searched for as functions of composition, energy, time, 
etc. Extensions to uncertainty in direction are trivial to include: simply divide a 
given event into a large number, M, say of subevents distributed appropriately and 
count each in the likelihood with a weight 1/M - this numerically convolves this 
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uncertainty with the hkehhood, and the size of any additional errors introduced by 
the procedure can be studied by varying M. 

5. Irresponsible Speculations 

It is interesting to speculate as to what any possible anisotropies might mean, so 
focussing on the idea of a few bright sources turning up, let us indulge in one wild 
idea, more for amusement than anything else! 

History has seen many attempts by human beings to either detect messages 
from extraterrestrial civilizations, or to send messages that might be received by 
intelligent life elsewhere. Carl Sagan, Linda Salzman Sagan, and Frank Drake^ 
had already proposed in 1972 that Pioneer 10, the first manmade craft to leave our 
solar system, should carry a message-bearing plaque in case some alien stumbled 
across it - and it did! This was a big step up from earlier proposals including 
one attributed to Carl Friedrich Gauss in the 1820's which suggested laying out 
geometrical patterns in the vast forests of Siberia by cutting our huge swaths of trees 
and planting wheat in their place. Austrian astronomer Johann Joseph von Littrow 
is alleged to have suggested the similarly non-environmentally friendly alternative 
of digging a vast circular canal (miles across) into the Sahara Desert, filling it with 
kerosene, and setting it on fire in the hopes that someone or something would notice 
and be impressed. 

Suggestions of radio-based communication go back to Nikola Tesla in the late 
1800's, and progressed through Frank Drake's 1960 "Project Ozma" to listen for 
signals at the 21 cm line of hydrogen to the modern SET I program. The interested 
reader can find a lot of historical data at the website E21 

Personally, I think it's interesting to wonder what would constitute an "impres- 
sive" signal to, or from, another civilization. By way of comparison, I am quite 
aware of the fact that gorillas have some degree of linguistic ability, but I am quite 
unlikely to spend much time trying having a sophisticated conversation with a go- 
rilla. While a gorilla and his friends may be quite impressed with themselves, the 
simple fact is that I'm unlikely to take much notice of their "intellectual" achieve- 
ments, and will think of them simply as apes. (There is, of course, no animosity 
intended - I'm just saying...) 

On the other hand, if a gorilla were to construct a laser pointer and fiash that 
around a bit, I might be a lot more impressed - at least here's a gorilla that has 
managed to make a controlled source of coherent light! In fact, the optical version of 
SETl^ makes the move away from looking for (or sending) signals by radio to using 
light. This, of course, offers the tremendous advantage that one can direct a laser 
towards a potentially inhabited planet and use a tiny fraction of the energy that 
would be required for a non-coherent source. Presumably that would be impressive 
to an alien civilization, and then they'd really want to talk to us - or would they? 

It's easy to argue that more exotic physics would be required for us to reach 
the threshold of being "interesting" . Perhaps we are awash in low energy neutrino 
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signals or gravitational waves for which we have no suitable detection devices. As far 
as we know, these might be possible channels for communication given a suitable 
trick of engineering - maybe we just haven't come up with it yet. It's certainly 
possible that the technologies needed to produce and detect these signals would 
develop hand-in-hand, in which case these channels would only be useful between 
civilizations at comparable levels of technology, and we would be able to participate 
in interstellar communications about the same time as we were able to listen in. 

UHECR's present a rather different kind of possible message-carrier: within our 
current knowledge of physics there's no way to make these with a device that would 
fit on our planet, while on the other hand, we are quite able to detect them! 

If it turns out that there is no natural mechanism of any kind to make UHECR's, 
and they must be created by artificial devices of a rather advanced technological 
nature, they would make great conversation openers from aliens who want to get 
our attention, even though we may not have comparable technology to signal back. 
Perhaps an anisotropy which revealed itself as pointing back to a planetary system 
could be a message to the effect of "Hey! Look at me! I've got advanced technology 
- advanced enough to make protons at 10^° eV and get them to you! These are 
pretty easy for you to detect (even though they're hard for me to produce), but 
at least you know I'm here!" In this sense, one might imagine that intelligently 
generated UHCER's could be cosmic signal flares to get our attention! 

This isn't really a completely serious suggestion, but then it's not entirely not 
serious. The only time it seems to have been discussed semi-seriously is in bars in 
Rio after this talk! 

6. Concluding Remarks 

I have reported on the first full-sky anisotropy search using data from the SUGAR 
and AGASA experiments. At present, low statistics and poor angular resolution 
limits our ability to perform a very sensitive survey, but we can at least have a 
preliminary look at the first moments in the angular power spectrum. The data 
are consistent with isotropy, though there appears to be a small excess for £ = 3, 
arising mostly from the SUGAR data. 

There are two caveats in this analysis which should be kept in mind. First, from 
the published SUGAR results, it is difficult to make an exact determination of the 
exposure, as the sensitive area of the experiment varied as a function of time. Here, 
we assumed an area-time product of approximately 775 km^ yr. Second, there is 
some uncertainty in the energy calibration. The SUGAR results are reported in 
terms of the number of vertical equivalent muons together with two possible models 
to convert this to primary energy. We have chosen the model yielding an energy 
spectrum which is in better agreement with the AGASA results It should be 
noted, though, that this spectrum does not agree well with the results of the High 
Resolution Fly's Eye (HiRes) experiment!^ Though there are uncertainties in the 
energy scale, the impact on this anisotropy search may not be so severe. This is 



5, 2008 6:33 WSPC/Trim Size: 9.75in x 6.5in for Proceedings 



main 



17 

because the energy cut of 10^^'^ eV is well above the last break in the spectral index 
at ~ lO"'^*'^ — 10^^ eV, and one would expect that all cosmic rays above this break 
share similar origins. 

In the n ear future we expect dramatically superior results from the Pierre Auger 
ObservatorjEll. This observatory is designed to measure the energy and arrival 
direction of ultra high energy cosmic rays with unprecedented precision. It will 
consist of two sites, one in the Northern hemisphere and one in the Southern, each 
covering an area S ~ 3000 km^. The Southern site is currently under construction 
while the Northern site is pending. Prelimi nar y results from the Southern site's 
engineering array are available in reference EHl Once complete, these two sites 
together will provide the full sky coverage and well matched exposures which are 
crucial for anisotropy analyses. Using Eq.|21 the calculated relative exposure is seen 
in Fig. 13 The base-line design of the detector includes a ground array consisting 




-100 -50 50 100 

Declination (deg.) 



Figure 5. Relative exposure of the full Pierre Auger Observatory. 

of 1600 water Cerenkov detectors overlooked by 4 fluorescence eyes. The angular 
and energy resolutions of the ground arrays are typically less than 1.5° (multi-pole 
expansion £ ~ 60) and less than 20%, respectively. The detectors are designed to be 
fully efficient {p « 1) out to ^^lnax — 60° beyond 10^^ eV, yielding a nearly uniform 
sky A w 1.4 X 10"* km^ sr.ESIln 10 yr of running the two arrays will collect w 4000 
events above Ei = 10^^'^ eV. As can be seen from Fig. 0] such statistics will allow 
one to discern asymmetries at the level of about 1 in 10^. 
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